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Summary. — We compare the luminosity function and rate inferred from the 
BATSE short hard bursts (SHBs) peak flux distribution with the redshift and lu- 
minosity distributions of SHBs observed by Smft/RETE II. The Swift/RETE II 
SHB sample is incompatible with SHB population that follows the star formation 
rate. However, it is compatible with a distribution of delay times after the SFR. 
This would be the case if SHBs are associated with the mergers of double neutron 
star (DNS) systems. DNS may be "primordial" or can form dynamically by binary 
exchange interaction in globular clusters during core collapse. The implied SHB 
rates that we find range from ~ 8 to ~ 30/i7 Gpc _3 yr _1 . This rate is a much higher 
than what was previously estimated and, when beaming is taken into account, it is 
comparable to the rate of neutron star mergers estimated from statistics of binary 
pulsars. If GRBs are produced in mergers the implied rate practically guarantees 
detection by LIGO II and possibly even by LIGO I. 

PACS 95.55.Ka - cosmology-observations. 
PACS 98.70.Rz - 7-ray sources; 7-ray bursts. 



1. Introduction 

In 1993 Koveliotou |l] noticed that the GRB distribution can be divided into two 
subsets of long and short bursts with a dividing duration of 2sec. As short bursts are 
also harder than long ones they are denoted as Short Hard Bursts (SHBs). 

The luminosity function (LF) and rate are fundamental quantities to understand the 
nature of these objects. Guetta & Piran (2005) (GP05) 4 L have shown how these quan- 
tities can be constrained by the observed peak flux distribution which is a convolution 
of these two unknown functions. 

Afterglows of SHB have only recently been observed by Swift |SJ E] and this has 
lead to the first identification of SHB host galaxies and redshift determination. In table 
1 we report the redshifts and luminosities of these bursts. As we can see from this table, 
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the average redshift of the short bursts is z ~ 0.4 much smaller than the average z of the 
long GRBs detected by Swift z ~ 2.3. 

As BATSE is less sensitive to short bursts than to long ones even an intrinsic 
SHB distribution that follows the star formation rate (SFR) gives rise to an observed 
distribution that is nearer to us. However the different sensitivity is not enough to 
explain the different z-distribution (z-DF) between long and short GRBs UJ. Therefore 
it is possible that SHBs do not follow SFR as the long bursts and that they are linked 
to the merger of double neutron star systems (DNS) ^Jj. In this case the SHB rate is 
given by the convolution of the star formation rate with the distribution P m (r) of the 
merging time delays r of the binary system. These delays reflect the time it takes to the 
system to merge due to emission of gravitational radiation. 

A delayed SFR distribution (that is intrinsically nearer) gives rise to an observed 
distribution that is nearer to us and that can explain the data for some choices of P(t). 
We consider two formation mechanisms. If two massive stars are born as a binary system, 
and the system remains bound after the super nova (SN) explosion of both components, 
a DNS is formed ("primordial" DNS). Another possibility is that at the moment of star 
formation the neutron stars (NSs) are not in the same binary system, but one of the NSs 
is in a binary with a low mass main sequence (MS) star. In globular clusters (GCs), such 
binaries are likely to have an exchange interaction with a single neutron star (see 
and references within). A significant fraction (~ 30%) of all NS mergers in the Universe 
may stem from such dynamically formed systems 

We show that the predicted z-DF is different for the two models, so that future 
z-observations may determine which formation channel (if any) is dominant. 



Table I. - The Swift/HETE II current sample of SHBs with a known redshift. 
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2. The luminosity function of the BATSE SHB sample 

Our data set and methodology follow |121 [S] . We consider all the SHBs detected 
while the BATSE onboard trigger ^2] was set for 5.5a over background in at least two 
detectors in the energy range 50-300keV. These constitute a group of 194 bursts. We 
assume the functional form of the rate of bursts (but not the amplitude) . We then search 
for a best fit LF. Using this LF we calculate the expected distribution of observed redshifts 
and we compare it with the present data. We consider the following cosmological rates: 

• (i) A rate that follows the SFR (We do not expect that this reflects the rate of 
SHBs but we include this case for comparison.). 

• (ii) A rate that follows the "primordial" DNs merger rate. This rate depends on 
the formation rate of NS binaries, that one can safely assume follows the SFR, 
and on the distribution of merging time delays, P m (r). This, in turn, depends 
on the distribution of initial orbital separation a between the two stars (r oc a 4 ) 
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and on the distribution of initial eccentricities. Both are unknown. From the co- 
alescence time distribution of six double neutron star binaries |14| it seems that 
P TO (log(r))dlog(r) ~ const, implying P m {r) oc l/r.[T5). Therefore, our best guess 
scenario is a SBH rate that follows the SFR with a logarithmic time delay distri- 
bution. 

• (iii) A rate that follows the SFR with a delay distribution P m (T)dT ~ const. 

• (iv) A constant rate (which is independent of rcdshift). 

For the SFR we employ Rs F2 of Porciani & Madau ^Sj : In models (ii) and (iii) the 
rate of SHBs is given by: 

/■*(*) 

(1) R SHB (Z)(X / dTR SF2 (t-T)P m { T ). 

Jo 



Table II. - Table 2: Best fit parameters Rate(z=0) , L* , a and f3 and their la confidence levels 
for models (i)-(iv). Also shown are the KS probability (pks) that the five bursts with a known 
redshift arise from this distribution. We show two results for KS tests one with GRB 0508132 
at z = 0.7 and the other at z — 1.8. Case ii a corresponds to case ii with an L* value lower by 
la than the best fit one. Other parameters have been best fitted for this fixed number. 
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Following Schmidt ^21 (see also we consider a broken power law peak LF with 

lower and upper limits, 1/Ai and A2, with power law indeces, a, f3 and luminosity break 
L*. 

We use Ai >2 = (30,100) 4 . In 9^ (denoted GP06 hereafter) we show that both 
limits are chosen in such a way that a very small fraction (less than 1%) of the observed 
bursts are outside these range. Hence one cannot infer anything from the observations 
on the LF in this range. Comparing the predicted distribution with the one observed by 
BATSE we obtain, the best fit parameters of each model and their standard deviation. 
The results are shown in table 2 and in figs 1 and 2 of GP06. 

3. — A Comparison with the current Swift-HET'E II SHB sample 

We can derive now the expected z-DF of the observed bursts' population in the 
different models. 

We assume that the minimal peak flux for detection for Swift is ~ 1 ph/cm 2 /sec 
like BATSE (note the different spectral windows of both detectors which makes Swift 
relatively less sensitive to short bursts). 

FigG] (left panel) depicts the expected observed integrated z-DF of SHBs in the dif- 
ferent models. As expected, a distribution that follows the SFR, (i), is ruled out by a KS 
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test with the current five bursts (pks < 1%)- This is not surprising as other indications, 
such as the association of some SHBs with elliptical galaxies suggest that SHBs are not 
associated with young stellar populations. 

A distribution that follows the SFR with a constant logarithmic delay distribution, 
(ii), is marginally consistent with the data (pks ~ 5 — 6%). The observed bursts are 
nearer (lower redshift) than expected from this distribution. If we use the Rowan- 
Robinson SFR 4 , rather than SF2 of Porciani-Madau, the situation is even more promis- 
ing pks ~ 10% (for either z=0.7 or z-1.8). When we move to a distribution that is ler 
away from the best fit distribution we find pks ~ 22 — 25% and even higher for the RR 
SFR. Thus the suggestion JHj that the observed sample rules out the NS merger model 
(with a logarithmic merger time distribution) was somewhat premature. Note however, 
that the local rate with this model (iv ) is sixteen times larger than the rate of the best 
fit model, (ii). This reflects the large flexibility in modeling the peak flux distribution. 

To demonstrate the flexibility of the data we have considered two other time delay 
distributions. A uniform time delay distribution (iii) and an overall constant SHB rate 
(iv) . Both models are compatible with the B ATSE SHB distribution and with the sample 
of SHBs with a known redshift (pxs ~ 80% and 40% respectively). This result is not 
surprising. The BATSE peak flux distribution depends on two unknown functions, the 
rate and the LF. There is enough freedom to chose one function (the rate) and fit for the 
other. 

In all models compatible with the five bursts with a known redshift, the intrinsic 
SHB rate is pushed towards lower redshifts. The inferred present rates, ~ 30, ~ 8 and 
~ 10/i7 Gpc _3 yr _1 for models (iii), (iv) and (ii CT ) respectively, are larger by a factor ten to 
fifty than those estimates earlier (assuming that SBHs follow the SFR with a logarithmic 
delay with the best fit parameters GP05). The corresponding "typical" luminosities, L*, 
ranges from 0.1 to 0.7 x 10 51 erg/sec. 

Let us consider now the possibility that short bursts come from dynamically formed 
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Fig. 1. - Left panel: A comparison between the expected integrated observed z-DF of SHBs 
for models (i)-(iv) and (iio-) and the distribution of known redshifts of SHBs. Right panel: A 
comparison between the expected integrated observed z-DF of SHBs for dynamically formed 
DNS 
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DNS [T§|. In this case, the delay time t is a sum of the time t cc until the dynamical 
formation of a DNS during core-collapse, and the time icw until the DNS merges. The 
resulting delay function (dp/dT)d yn of dynamically formed DNSs is given by 



(2) (f\ =± fdtj^{t cc) r cc dt GW d J^ {tGW ). 

\ dT / dyn dT JO *cc Jo *GW 

The resulting delay function is shown in Fig. 2 of [PQ. Fig.l (right panel) shows the 
comparison between the expected and observed integrates z-DF. Our model for dynamical 
mergers fits the observed SHBs distribution much better (pk s ~ 70% — 80%) than the 
model for primordial distribution with P(r) ~ 1/r (pks ~ 5% — 20%) . 



4. Conclusions and Implications 

We have repeated the analysis of fitting the BATSE SHB data to a model of the 
luminosities and rates distributions. Our best fit logarithmic distribution model is similar 
to the best fit logarithmic model presented in GP05. A main new ingredient of this work 
is the fact that we consider several other models. We confirm our earlier finding that the 
BATSE data allows a lot of flexibility in the combination of the rates and luminosities. 

A second new ingredient of this work is the comparison of the best fit models to 
the small sample of five Sm/i/HETE II SHBs. The Swift/EETE II data gives a new 
constraint. This constraint favors a population of SHBs with a lower intrinsic luminosity 
and hence a nearer observed z-DF. It implies a significantly higher local SHB rate - a 
factor of ten to fifty higher than earlier estimates. The new observations of Swift show 
that the SHBs are nearer than what was expected before and therefore, their luminosity 
is lower and their local rate is higher. We stress that this new result was within the la 
error of the model presented in GP05, which had a very large range of allowed local rates 
and typical luminosities. 

Provided that the basic model is correct and we are not mislead by statistical (small 
numbers), observational (selection effects and threshold estimates) of intrinsic (two SHB 
population) factors we can proceed and compare the inferred SHB rate with the obser- 
vationally inferred rate of NS-NS mergers in our galaxy [201 ■ This rate was recently 
reevaluated with the discovery of PSR J0373-3039 to be rather large as 80lgg°/Myr. 
Although the estimate contains a fair amount of uncertainty 2b . If we assume that this 
rate is typical and that the number density of galaxies is ~ 10~ 2 /Mpc 3 , we find a merger 
rate of SOOt^/GpcVyr. Using a beaming factor of 30-50 for short bursts [7] this rate 
implies a total merger rate of ~ 240 — 1500/Gpc 3 /yr for the three cases (iii), (iv) and 
(iier). The agreement between the completely different estimates is surprising and could 
be completely coincidental as both estimates are based on very few events. 

If correct these estimates are excellent news for gravitational radiation searches, for 
which neutron star mergers are prime targets. They imply that the recently updates 
high merger rate, that depends mostly on one object, PSR J0737+3039, is valid. These 
estimate implies one merger event within ~ 70Mpc per year and one merger accompanied 
with a SHB within ~ 230Mpc. These ranges are almost within the capability of LIGO I 
and certainly within the capability of LIGO II. If correct these estimates of the rate are 
excellent news for gravitational radiation searches, for which neutron star mergers are 
prime targets. 
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We have considered also the possibility that the DN system is formed dynamically 
and we have shown that the z-DF of short GRBs expected from dynamical formation 
and subsequent merger of DNSs is markedly different from that expected from primordial 
DNS mergers (under the assumption that the distribution of delay times for primordial 
DNSs follows (dp/dr) P ii m l/ r )- The large time for core collapse shifts the distribution 
of dynamically formed DNS mergers to low redshift. The observed z-DF of short GRBs is 
consistent with that expected for mergers of dynamically formed DNSs, and appears to be 
inconsistent with that expected for primordial DNS mergers if the best fit parameters are 
considered We note, however, that current data do not allow to rule out a z-DF consistent 
with that expected for primordial DNS mergers, since redshifts were obtained only for 
a minority of the detected short GRBs. This may be due to a bias against obtaining 
redshift information for high redshift (faint) short GRBs (see Fig. 4, table 1 and discussion 
at the end of §3.3 of ^Hl). Future observations should allow to better constrain the z-DF 
of short GRBs, and thus to differentiate between models. For example, detection of only 
a few high redshift (z > 2) short GRBs would severely constrain the contribution of 
dynamically formed DNSs. 

To conclude we stress that we have assumed that the luminosity function has a lower 
limit of L*/30. This was just because even if such a limit does not exists weaker bursts 
would be barely detected. The current peak flux distribution of BATSE burst cannot 
confirm (or rule out) the existence of such population (note however, that Tanvir et al., 
|22| suggest that such a population exists on the basis of the angular distribution of 
BATSE SHBs). If such weak bursts exist then, of course, the overall merger rate will be 
much larger Such events will provide such a high rate that soon LIGO I will begin 
to constrain this possibility. 
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